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Generating dynamical black hole solutions 
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We prove a theorem that characterizes a large family of non-static solutions to Einstein equations, 
representing, in general, spherically symmetric Type II fluid. It is shown that the best known 
dynamical black hole solutions to Einstein equations are particular cases from this family. Thus we 
extend a recent work of Salgado [I^] to non-static case. The spherically symmetric static black hole 
solutions, for Type I fluid, are also retrieved. 
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I. INTRODUCTION 



One of the most famous conjectures in general relativity is the so called cosmic censorship conjecture (CCC) [l| 
which states that, for physically reasonable initial data, space-time cannot evolve towards a naked singularity, i.e., the 
space-time singularity is always hidden inside black holes, indicating that a far away observer will not be influenced 
by it. Despite almost 30 years of effort we are far from a general proof of the CCC. In fact, it turns out that such 
a theorem is intractable due to complexity of the Einstein field equations, hence metrics with special symmetries 
— i' are used to construct gravitational collapse solutions. As a result, there are very few inhomogeneous and nonstatic 
^ , solutions known, one of them is the Vaidya metric. 
The Vaidya metric [21] , which has the form 



ds^ 



2m{v) 



dv^ + 2edvdr + r^{d9^ + sin^ed<l)^), e±l (1) 



^ ' is a solution of Einstein's equations with spherical symmetry for a null fluid (radiation) source described by energy 

^^ , momentum tensor Tab — ip^ah, la being a null vector field. For the case of an ingoing radial flow, e = 1 and m{v) is 

■^ ■ a monotone increasing mass function in the advanced time v, while e — —1 corresponds to an outgoing radial flow, 

_. -j with m(v) being in this case a monotone decreasing mass function in the retarded time v. The Vaidya's radiating 

star metric is today commonly used for two purposes: (i) As a testing ground for various formulations of the CCC. 

f^ (ii) As an exterior solution for models of objects consisting of heat-conducting matter. Recently, it has also proved to 

be useful in the study of Hawking radiation, the process of black- hole evaporation ^ , and in the stochastic gravity 

program [J]. Also, several solutions in which the source is a mixture of a perfect fluid and null radiation have been 

obtained in later years [3|. This includes the Bonnor- Vaidya solution la| for the charge case, the Husain solution ig| 

with an equation of state P = kp. Glass and Krisch [8] further generalized the Vaidya solution to include a string 

fluid, while charged strange quark fluid (SQM) together with the Vaidya null radiation has been obtained by Harko 

^ and Cheng [9|]. Wang and Wu |lfl] further extrapolated the Vaidya solution to more general case, which include a 

C^ [ large family of known solutions. Their generalization comes from the fact that the linear superposition of particular 

solutions is also a solution of Einstein's equations. The Vaidya solution in the brane-world was discovered by Dadhich 

and Ghosh [IH 

Recently, Salgado |12l | proved a theorem characterizing a three parameter family of static and spherically symmetric 
solutions (black holes) to Einstein equations by imposing certain conditions on the energy momentum tensor (EMT) 
(see also [13nl5]). His work was extended to higher dimensional spacetime by Gallo [16]. However this is obviously 
not the most physical scenario and one would like to consider dynamical black hole solutions, i.e., black holes with 
non-trivial time dependence. 

In this paper, we consider an extension of Salgado's work, so that a large family of exact spherically symmetric 
dynamical black hole solutions are possible, including it's generalization to asymptotically de Sitter/ Anti-de Sitter 
(dS/AdS). It turns out that the family of solutions generated represent generalization of Vaidya solutions that include 
most of the known solutions. 
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II. THE DYNAMICAL BLACK-HOLE SOLUTIONS 



Theorem - I; Let (M, gat) be a four dimensional space-time [sign{gab) = (—,+,+,+)/ such that (i) It is non-static 
and spherically symmetric, (ii) it satisfies Einstein field equations, (Hi) in the Eddington-Bondi coordinates where 
ds^ = —A{v,r)'^f(v,r) dv^ + 2eA(w,r) dv dr + r^idO"^ + sin'^OdcjP), the energy-momentum tensor T""^ satisfies the 
conditions T^ = T^ and Tg = kT^ , {k = const, e TZ) (iv) it possesses a regular Killing horizon or a regular origin. 
Then the metric of the space-time is given by 



where 



ds' ^- 



1- 



2m{v, 



m(v, r) 



dv^ + 2edvdr + r^{d0^ + sin^Odcj}^), 



(e = ±l) 



' M{v 

^ M{v) -47rC(w)lnr 



47rC(-») ^2fc+l 



if C{v) = 0, 

if C{v) 7^ and /c =^ -1/2, 

if C{v) ^ and /c = -1/2. 



(2) 



(3) 



Tb=-S^dmg[l,l,k,k], 



and 



rjir . 



1 dM 

TttF^ dv 

1 dM 

47rr^ dv 



1 ac 2fe-i 

2fe+l dv ' 


iik^ 


-1/2, 


y^^-r 


if fc = 


-1/2. 



(4) 



(5) 



Here, M(v) and C(v) are the integration constants whose values depend on the boundary conditions and the 
fundamental constants of the underlying matter. 

Proof: Expressed in terms of Eddington coordinate, the metric of general spherically symmetric space-time [17| is, 

ds^ = -A{v, rff{v, r) dv^ + 2eA{v, r) dv dr + r^{d9^ + sin^edtj)^). (6) 



Here A{v,r) is an arbitrary function. It is useful to introduce a local mass function m{v,r) defined by f{v,r) = 
1 — 2m{v,r)/r. For m{v,r) = m{v) and A — 1, the metric reduces to the standard Vaidya metric. It is the field 
equation G^ = that leads to A{v,r) = g{v). However, by introducing another null coordinate TJ = J g{v)dv, we can 
always set without the loss of generality, A{v, r) = I. Hence, the metric takes the form. 



ds^ = - 



1- 



2m{v, r) 



dv^ -ir2edvdr + r^{de'^ 



'r 



(7) 



Therefore the entire family of solutions we are searching for is determined by a single function m{v,r). Henceforth, 
we adopt here a method similar to Salgado [I4I which we modify here to accommodate the non static case. In what 
follows, we shall consider e = 1. The non- vanishing components of the Einstein tensor are 

(8a) 

(8b) 

(8c) 



where {a;"} = {u, r, 9, (/p}. The Einstein field equations are 



Ra 



Gl 


2 dm 
r^ dv ' 




g: 


= g; = - 


2 dm 

J.2 Qj. 


Gl 
uat 


= C^ = - 
ons are 


Id^m 
r dr"^ 


ab - 


- -^RQab = 


SnTab , 



(9) 



and combining Eqs. ([5]) and ©, we have if a ^ &, T^° = except for a non-zero off-diagonal component T^. It may 
be recalled that EMT of a Type II fiuid has a double null eigen vector, whereas an EMT of a Type I fiuid has only 



one time-like eigen vector '2]^. In addition, we observe that the metric ([7]) requires that T^ = T^. Thus the EMT 
can be written as : 



^b — 



fT^ \ 

t; t; 

T^ 

V T^ J 



Enforcing the conservation laws '^ aT^ — 0, yields the following non-trivial differential equations: 



dr 


r 


-T'el 


dv 


dr 


2 



Using the hypothesis that Tg — kT^, we obtain the following linear differential equation 



or r 



which can be easily integrated to give 



rpr 



C{v) 



j,2(l~k) ' 

where C{v) is an integration constant. Then, using hypothesis (iii), we conclude that 

Civ) 



rpa 

f" ~ j,2(l-k) 



diag[l, l,fc, k]. 



Now using Eqs. (|8]) and ([9|) and ([13]), we get dm/dr — —4:nC{v)/r "^^ , which trivially integrates to 

f M{v) if C{v) = 0, 

M{v) - i|g^r2'=+i if C{v) ^Oemdk^ -1/2, 

^ M{v) - 4:nC{v) Inr if C{v) ^ and k = -1/2. 



m(v, r) — < 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 



Here the function AI{v) arises as an integration constant. What remains to be calculated is the only non-zero off- 
diagonal component T^ of the EMT. From Eqs. ([8]) and ([9]), one gets 



j.r 



1 dm 

Airr^ dv ' 



which, on using Eq. (llSp . gives 



rpr 



1 8M 
47rr^ dv 


1 dC 2k- 
2k+l dv ' 


'I 


1 dM 
47rr^ dv 


-;^^l- 





if fc ^ -1/2, 
if fc = -1/2. 



(16) 



(17) 



It is seen that Eq. pT|) is identically satisfied. Hence the theorem is proved. 

The theorem proved above represents a general class of non-static, spherically symmetric solutions to Einstein's 
equations describing radiating black-holes with the EMT, which satisfies the conditions in accordance with hypothesis 
(iii). The static solutions derived in [l^] can be recovered by setting M{v) = M, C{v) — C, with M and C as 
constants, in which case matter is Type I. The solutions generated here highly rely on the assumption (iii). On the 
other hand, although hypothesis (iv) is not used a priori for pro ving the result, but it is indeed suggested by regularity 
of the solution at the origin, from which, T^ = T^\r=o (see [13| for further details). 



TABLE I: Energy-momentum tensor and the corresponding non-static space-time they generate alongwith the functions M{v) 
and C{v), and the fc-index associated with each space-time. 



Energy-momentum tensor 



Space-time 



Functions: A/(u)and C{v) 



fc-index 



rj-^a n. rpr 1 d hi 



rpV rpr ^ a 



n = -|^ diag[l, 1, -1, -1] 

rpr _ 1 r rdM 



1^ 



Ta 
h 



g(f) 



.... T^T+TT diag[l, 1, -m, -m] 
rpr _ 1 \ai 1 ag] 



Vaidya 

dS/adS 

Global monopole 

Bonnor- Vaidya 
Husian 



M{v), C{v) = 
M{v) = 0, C{v) = -A 
M{v) = 0, C{v) = -^ 






fc = 


1 


fc = 





fc = 


-1 


k = 


-m 



kT, 



(f)r- 



Then the theorem I 



The family of the solutions outlined here contains, for instance, Bonnor- Vaidya, dS/AdS [To|, global monopole jTsI, 
Husain, Harko-Cheng SQM solution [3,|l3: and Glass-Krisch string solutions [1, [20l- Obviously, by proper choice of 
the functions M{v) and C[v), and A;— index, one can generate as many solutions as required. The above solutions 
include most of the known spherically symmetric solutions of the Einstein field equations. Some of the examples of 
EMT satisfying the conditions of the theorem and which generates the known space-times are summarized in Table L 
The solutions discussed in the section are characterized by two arbitrary functions M(v) and C(w), and the cosmo- 
logical constant A. Thus one would like to generalize the above theorem to include A. We can show that the energy 
momentum tensor components, in general, can be written as, T^ = T^f\f, — ^S^ 0, 13, where A is the cosmological 

constant and T^rs^^ is energy momentum tensor of the matter fields that satisfy T^^.g 
can be generalized as: 

Theorem - II; Let (M, gab) be a four dimensional space-time [sign{gab) = (^, +, +7 +)/ such that (i) It is non-static 
and spherically symmetric, (ii) it satisfies Einstein field equations, (Hi) the total energy-momentum tensor is given by 
T^ = T?'r-,^ — -^^b' '"^here A is the cosmological constant and T?'rs^ is energy momentum tensor of the matter fields, (iv) 
in the Eddington coordinates where ds^ = —A{v, r)'^f{v, r) dv"^ -\- 2A{v, r) dv dr -\- r'^{d9'^ -\- sin^Odcj)^), the EMT T?-r-,^ 
satisfies the conditions Th^^ — TL^^, T?r^g — kTL^^, [k — const. £ TZ), (v) it possesses a regular Killing horizon or 
a regular origin. Then the metric of the space-time is given by metric (0), where 

r M{v) + ^ if C{v) = 0, 



m{v, r) 



M{v) 



47rC(t>) 2fc+l 
2k+l ' 



. M{v) -47rC(w)lnr-t- 



6 

Arf. 
6 



if C{v) y^ and fc =^ -1/2, 
if C{v) ^ and fc == -1/2. 



(18) 



and 



TS = ;|(Sydiag[l, 1, fc, k] - ^cliag[l, 1, 1, 1] 



y = 



1 dM 

TttF^ dv 



1 dC 2k~l 
2fe+l dv ' 



1 dM _ I dC_Y-^^ 
A-KT-^ dv r^ dv 



iik^ -1/2, 
if fc== -1/2. 



(19) 



(20) 



Here, M(v) and C(v) are the integration constants whose values depend on the boundary conditions and the 
fundamental constants of the underlying matter. 



Again the conservation law V aT^ — leads to 



dT, 



U)r 



dr 



'in 



f)r 



^(/)e)- 



(21) 



Using the assumption made above that Tf^^g = kTLs, we obtain 

^---(l-fc)^(7).' (22) 

which can be easily integrated to give 

_ Cjv) 

where C{v) is an integration constant. Then, using hypothesis (Hi) and (iv), we conclude that 

^b = ^S^diagil, 1, fc, k] - ^diag[l, 1, 1, 1]. (24) 

Now using Eqs. dH), ^ and ([Ml), we get 



dm 
dr 



= -47rr^ 



Civ) A 



(25) 



j.2(l-k) g^ 

which, on integration, gives, (TTSl) . Next, to calcualte the component T^ of the EMT. From Eqs. ([8]) and ^, one gets, 
which, on using Eq. ([T8|) . gives ([20| . Hence the theorem is proved. 



III. ENERGY CONDITIONS 

The family of solutions discussed here, in general, belon gs t o Type II fluid defined in [2l|. When m = m{r), we 
have /x=0, and the matter field degenerates to type I fluid jlO|. In the rest frame associated with the observer, the 
energy-density of the matter will be given by (assuming A = 0), 

Ai = r;, p = -t,* = -t;^--5^^3^, (27) 

and the principal pressures are Pi = T!* (no sum convention). Therefore Pr = T^ = —p and Pe ^ Pip = kPr = —kp 
(hypothesis {in)). 

a) The weak energy conditions (WEC): The energy momentum tensor obeys inequality Tab'w°"w^ > for any timelike 
vector, i.e., 

M > 0, p > 0, Pe > 0, Py > 0. (28) 

We say that strong energy condition (SEC), holds for Type II fluid if, Eq. (1^ is true., i.e., both WEC and SEC, for 
a Type II fluid, are identical. 

b) The dominant energy conditions : For any timelike vector Wa, T°'^'WaWb > 0, and T'^^Wa is non-spacelike vector, 
i.e., 

M>0, p>Pe,P^>0. (29) 

Clearly, (a) is satisfied if C{v) < 0, fc < 0. However, /i > gives the restriction on the choice of the functions M{v) 
and C{v). From Eq. ([5]), (fc ^ —1/2), we observe p > requires, 

dM 47r dC o,,, , „ 

-d^~WTT)^^">'- (^«) 

This, in general, is satisfied, if 

dM dC dC 

-—- > 0, and, either -- > and fc < -1/2, or ^— < and fc > -1/2. (31) 

ov ov ov 

On the other hand, for fc = -1/2, /i > if dM/dC > Anlnr. The DEC holds if C(v) < and -1 < fc < 0, and the 
function M is subject to the condition (|31l) . Clearly, < — fc < 1. 



IV. SINGULARITY AND HORIZONS 



In this section, we shall discuss the physical properties of the solutions. The Ricci R = RabR""^ , Rab the Ricci tensor 
and Kretschniann invariants(K = RabcdR'^'^'^'^ , Rabcd the Riemann tensor), for the metric ^, reduces to: 



R = 



1287r2C(w)2(l + fc2) 327rAC(w)(l + /s) 



j.A(l~k) 



j.2{l-k) 



4A2. 



(32) 



K 



48M-^ 



2567r2C(w)= 



„2k+l 



2fc + l 



[(/c-l)2(2fc+l)2 + 2fc(4fc-l)] 



lQnM{v)C{vy^+l^^ _ ^^^ _ g^2] ^ 64^AC(«) 



r6(2fc+l) 
x[(l-4fc) + (l-4fc2)(fc + 2)] 



3r6 



^fc+2 



2fc + l 



if k^ -1/2, 



(33) 



K = 



48M2 



647rC(w) 



-47rC(w)[5-31nr]lnr 



137rC(w) + M(f)[5-61nr] 
8A3 32AC(u)7rr3 



(34) 



6 



if 



-1/2. 



These invariants are regular everywhere except at the origin r — Q, where they diverge. Hence, the space-time has the 
scalar polynomial singularity [2l| at r = 0. The nature (a naked singularity or a black hole) of the singularity can be 
characterized by the existence of radial null geodesies emerging from the singularity. The singularity is at least locally 
naked if there exist such geodesies, and if no such geodesies exist, it is a black hole. The study of causal structure of 
the space-time is beyond the scope of this paper and will be discussed elsewhere. 

In order to further discuss the physical nature of our solutions, we introduce their kinematical parameters. Following 
York [22| a null- vector decomposition of the metric ([2]) is made of the form 



gab 



where, 



8"" I 



nJb - Lm + "tab 

2m{v, r) 



51 + s: 



l"" lab ^ 0;iab n'' ^ 0, 



(35) 

(36a) 
(36b) 

(36c) 



with m{v,r) given by Eq. (|3]). The optical behavior of null geodesies congruences is governed by the Raychaudhuri 
equation 



dv 2 



-WafcCJ 



ab 



(37) 



with expansion 0, twist w, shear a, and surface gravity IC. The expansion of the null rays parameterized by v is given 
by 



where the V is the covariant derivative. In the present case, tr = w = [.22j |. and the surface gravity is, 

K. = -n'^l^Vbla- 



(38) 



(39) 



As demonstrated by York l22l|, horizons can be obtained by noting that (i) apparent horizons are defined as surface 
such that 8 ~ and (ii) event horizons are surfaces such that dQ/dv c± 0. Substituting Eqs. (|18p . ([551 and ([5^ into 
Eq. dSll), we get, (fc ^ -1/2) 



9^1 

r 



l-^^^+Q2(„)^2fe_^2^2 



(40) 



where x^ = A/3, and Q'^{v) = 8TrC{v)/2k + 1. Since the York conditions require that at apparent horizons 9 vanish, 
it follows form the Eq. (HUl) that apparent horizons will satisfy 

xV^ - gV^'^+i - r + 2M = 0, (41) 

which in general has two positive solutions. For X^ — Q^ = Oi '^6 have Schwarzschild horizon r — 2M, and for 
M = Q^ = we have de Sitter horizon r = 1/x- As mentioned above, for k = —1, one gets Bonnor-Vaidya solution, 
in which case the various horizons are identified and analyzed by Mallett [2^1 and hence, to conserve space, we shall 
avoid the repetition of same. For general k, as it stands, Eq. (|4T|) will not admit simple closed form solutions. However, 
for 



l2 ^2 1 



Q=Qi 



i2k + l) 



2k "'^'^ 



2M(2fc + l) 



with X — 0, the two roots of the Eq. (1411) coincide and there is only one horizon 

2M(2fc + 1) 



r = 



2fc 



(42) 



(43) 



For Q^ < Qc there are two horizons, namely a cosmological horizon and a black hole horizon. On the other hand if, 
the inequality is reversed, Q^ > Q\ no horizon would form. 

V. CONCLUDING REMARKS 

In conclusion, we have extended to non-static case a recent theorem [12] and it's trivial extension (that includes 
cosmological term A), which, with certain restrictions on the EMT, characterizes a large family of dynamical black 
hole solutions, representing, in general, spherically symmetric Type II fluid. The solutions depend on one parameter 
k, and two arbitrary functions M{v) and C{v) (modulo energy conditions). It is possible to generate various solutions 
by proper choice of these functions and parameter k. Many known solutions are identified as particular case of this 
family and hence there exists realistic matter that follows the restrictions of the theorem. 

The family of solutions discussed here, in general, belongs to Type II fluid. However, if M = C ~ constant, we 
have fi—O, and the matter field degenerates to type I fluid and we can generate static black hole solutions by proper 
choice of these constants. 

A rigorous formulation and proof for either version of CCC is not available. Hence, examples showing occurrence 
of naked singularities remain the only tool to study the various aspects of CCC. However, the lack of exact solutions 
of Einstein field equation makes it very difficult. As a consequence, we are far from complete understanding of CCC 
even in the simple case of spherical symmetry. The solutions presented here can be useful to get insights into more 
general gravitational collapse situations and in general better understanding of CCC that may help to put CCC in 
precise mathematical form. 
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